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We study the one-loop potential of large- N QCD with adjoint Dirac fermions. Space-time 
is a discretization of R 3 x S 1 where the compact direction consists of a single lattice site. We 
use Wilson fermions with different values of the quark mass m and set the lattice spacings in 
the compact and non-compact directions to be a t and a s respectively. Extending the results of 
JHEP 0906:091,2009, we prove that if the ratio £ = a s /at obeys < £ < 2, then the minimum 
of the one-loop lattice potential for one or more Dirac flavors is symmetric at the chiral 
point. For £ = our formulas reduce to those obtained in a continuum regularization of the 
R 3 , and our proof holds in that case as well. As we increase m from zero, we find a cascade 
of transitions where Zn breaks to Zk- For very small masses, K ~ l/(atm) 3> 1, while for 
large masses K ~ 0(1). Despite certain UV sensitivities of the lattice one-loop potential, 
this phase structure is similar to the one obtained in the continuum works of Kovtun-Unsal- 
Yaffe, Myers- Ogilvie, and Hollowood-Myers. We explain the physical origin of the cascade of 
transitions and its relation to the embedding of space-time into color space. 



1 



I. INTRODUCTION 



There is much recent interest in the phenomenon of space-time reduction of SU(N) gauge 
theories [l). Reasons for this include the hope to understand nonperturbative dynamics of 
QCD with analytic small volume techniques [2], and the possibility to reduce the numerical 
cost in calculating properties of large- N gauge theories on the lattice. For related works on 
non-lattice simulations see relevant references in jsj. 

The nonperturbative definition of spacetime reduction requires a nonperturbative definition 
of the gauge theory whose volume one wishes to reduce. This requirement naturally leads 
one to a lattice definition. Indeed, all the known proofs of spacetime reduction that go 
beyond perturbation theory involve a lattice regularization. These proofs (see for example 



3na perti 

5. QUE 



Refs. Hl,y,l5|) state the following:. 

Given an SU (N) lattice gauge theory defined by a set of dimensionless couplings (the lattice 
't Hooft coupling g 2 N, a set of bare masses in lattice units ami, a>m-ii • • etc.), and residing on 
a geometry defined by a set of integers £1,2,... counting the number of sites along the different 
compact directions, then, under certain conditions, physical observables (e.g. Wilson loops, 
meson propagators, quark condensates, etc.) are independent of £1,2,... when N = 00. These 
conditions are 

1. The ground state is translation invariant. 

2. The ground state is invariant under center transformations that rotate Polyakov loops 
that wrap the compact directions. 

3. The ground state obeys cluster decomposition (large- N factorization). 

The validity condition that is hardest to fulfill seems to be no. (2). In the theory of 
interest to the current paper - four dimensional SU(N) gauge theory defined on R 3 x S 1 with 
adjoint fermions - this condition involves the center symmetry Zn whose order parameter is 
the Polyakov loop Q that wraps the S 1 . The ultimate check of whether the ground state of 
this theory respects the center symmetry can be performed only with nonperturbative lattice 
simulations. Nevertheless, one can approach this question analytically in certain regimes. 
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In this paper we focus on the regime of weak-coupling and use perturbation theory to 
calculate the effective potential that governs the expectation value of Q. This analysis is 
justified when the inverse length of the S 1 is large compared to Aqcd, and, for the massless 
case, the calculation was performed in the continuum theory by Kovtun, Unsal and Yaffe in 
Ref . ((J . For nonzero masses, the continuum calculation was done by Myers and Ogilvie, and 
Hollowood and Myers in Refs. jl, 8]. 1 

The picture emerging from these weak coupling calculations is interesting 2 : At zero mass 
the ground state is Zn invariant. As the mass increases from zero there appears a cascade 
of transitions where the center breaks to smaller and smaller subgroups. When the mass is 
infinite the symmetry breaks completely. 3 

The aim of this paper is two-fold. First, we wish to gain physical insight on the origin of 
this cascade of transitions. Second, we wish to generalize the calculations of Refs. js S] to the 
lattice regularization defined with a single site in the compact direction. Our interest in this 
setup derives from the UV sensitivities of the one-loop lattice potential that we analyzed in 



Ref. 



10] . These can make the Zn realization of the lattice theory different than that of the 



continuum theory, and even dependent on the lattice action. Therefore, we choose an action 
which is the one simulated in Ref. [lsJ - pure Yang-Mills fields regularized with the simple 
plaquette Wilson action and a fermion regularization of the Wilson type. We also allow for 



1 While Ref. [7[ preceded Ref. [8j , the latter used techniques that allow one to avoid a numerical minimization 
of the one-loop potential. Also, we note in passing that while some of the results of Ref. 8] involve the 
space S 3 x S 1 , then it is the result in the large S 3 limit, where the sphere becomes an R 3 , that is relevant 
to our calculation. For earlier results on S 3 x S 1 see Ref. Q. 

2 Here we restrict to more than a single Majorana fermion. The SUSY case is more complicated and we refer 
to Ref. [6| for further discussion. 

3 Clearly this assumes that N is non-prime (otherwise Zn would not accommodate subgroups). If N is prime, 
then a would-be instability towards a ground state that has a Zk symmetry with K < N, will actually lead 
to a ground state that completely breaks Zn, but that 'looks' approximately Zk symmetric. For example, 
suppose N is odd and choose the mass to be such that for an SU(N') gauge group with N' = N + 1 the 
center Zn> would break to Z 2 - This will make the eigenvalue density of the Polyakov loop double peaked. 
The intact Z 2 symmetry will ensure that the two peaks have identical heights. For the same mass the 
eigenvalue density in the SU (N) case will also have two peaks, but of different heights. At large enough 
values of N, however, the difference between the SU(N) and SU(N + 1) theories will be washed away 
and the peaks will differ by a 1/N amount. Thus, even the N = odd case can experience an approximate 
Zn —s- Z 2 symmetry breakdown. 
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an asymmetry between the lattice spacing in the non-compact directions, a s , and the lattice 
spacing a* in the compact direction. The ratio of the two spacings is denoted by 

£ = a s /a t . (1.1) 



The current paper can be viewed as a continuation of Ref. 



101 ] . where we analyzed the 



same setup studied here, but in a more restricted way and with somewhat different methods. 



In particular, in Ref. [10[ we compared the energy of different vacua by explicitly studying 
particular realizations of the eigenvalues of Q. In contrast, in the current paper we directly 
investigate the eigenvalue density of Q and the way it may realize different symmetry break- 



downs. This method is advantageous to the one used in Ref. [10| because it allows one to 
study many type of symmetry realizations. Indeed, with the relatively inefficient calculation 
method of Ref. jioj, we only examined the energy of the ground state where Z^ is intact, of 
the one that completely breaks Z^, and of the one that partially breaks Zn to Z2. In the 
current paper we study many more possible symmetry realizations of the form Zn — » Zk with 
K < N, and for simplicity assume that N/K = integer (but see footnote [3]) . 

Our current calculation method also allows us to make certain analytic statements that do 
not rely on numerical evaluations of the one-loop free energy, and leads to a simple picture of 
the emergent phase structure that is connected to the embedding of spacetime in color space. 
The following is the outline of this paper. In Section [TT] we show how to recast the one- 



loop potential V(Q) calculated in Ref. 10] in a form useful to the purposes of this paper. 
Specifically, we write 

00 

V(n) = 2j2V r I trfi r | 2 + constant. (1.2) 

r=l 

Since the density p of the eigenvalues of fl is defined by 

pb») = (jf'E HP -<*)), (1-3) 



a=l I 

then its r th moment, p r , defined by 

P, = Re J^e^p(u), (1.4) 

obeys 

p r = Re (^trfi r V (1.5) 



This means that the sign of V r determines whether the eigenvalue density develops an insta- 
bility in its r th moment and breaks the symmetry from Zn to Z r . 

In Section II III we minimize V(Q) by investigating the signs of V r for different values of 
m,a s ,a t , and Nf. The mathematical form we obtain for V r makes clear the spacetime em- 
bedding into color space, and leads one to anticipate and understand the transition cascade 

nfi 

obtained in Refs. [7|, |8|. We discuss this point in Section [II Bl Then, we turn to minimize the 
one-loop potential. We start with the m = case, where we determine the sign of V r analyti- 
cally: for finite values of a S) t see Section IIII Al and for the a s — > limit see Section IIIIBl We 
proceed in Section HlI CI to calculate V r numerically for general values of the parameters (since 
V r takes the form of an integral over complicated functions defined within the lattice spatial 
Brillouin zone, we cannot evaluate it analytically for generic parameters). We summarize our 
results in Section HVl 



We note in passing that in a related study, Ref. calculated the one-loop potential of the 
single site theory that we study here, and regulated the UV of the R 3 in the continuum. As 



Ref. lOf showed, if certain relevant operators are added to the one- loop potential of Ref. 11] , 
then such a continuum analysis is expected to yield results which are consistent with those 
we report here (at least up to finite lattice corrections) and with those reported in Ref. 6]. 



Nevertheless, the results of Ref. Uj indicated that the zero mass system has a ground state 
that breaks the Zjy symmetry; a result that contrasts the results of the works mentioned 
above. Our current analysis shows that this Z^ breaking is a result of an erroneous sign in 



the final expression for V(Q) that Ref. 11] minimized. For clarity, we re-derive in Appendix |D1 
the form of V(O) with a continuum regulator for the non-compact UV and show that it yields 
results which are consistent with those appearing in Ref. |6| and in the current paper. 

II. THE ONE LOOP POTENTIAL 

In this section we calculate the one-loop potential and recast it in a useful form for our 
purposes. This is done in Section Hi Al In Section Hi Bl we discuss some of its general features. 
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A. The calculation 



In Ref. [10[ we calculated the one-loop potential V as a function of the Polyakov loop Q 



that wraps the S 1 . In terms of the eigenvalues e l9a of Q the result was 



-2JV) log 



— p + 4 sin 



1 i2 



£ 2 



p +sin 2 (e a -e b )+m 2 w (p,9 a -e 



(2.1) 



To obtain this result we neglected the contribution of a set of fluctuations in the N diagonal 
components of the gauge fields, which is justified at large- N. 
The quantities p and p are the lattice momenta 



p 2 = Y 4 sin2 fe/ 2 ) > p = X! sin2 



(2.2) 



4=1 



and miy accounts for the contribution of the bare mass and the Wilson term to the propagator 



m w = at~m + — p 2 + 2 sin 2 



To proceed we write 



(2.3) 



with 



V(Q) = 2 Re E^e" 

aj^b r=l 



+ constant 



2 V) V r |trfi r | 2 + constant, 

r=l 



do; 



K = Re — e J ™ 



2tt 
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dp 



2tt 
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1 a2 



p + sin (w) + 77% (p, u) 



In Appendix |A] we show that V r is given by 

3 



1 



dp 



r J-7T V 27r 



2iVy e ~ rE ^) — e~ rEa ^ 



r > 1, 



(2.4) 



(2.5) 



(2.6) 
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with Ep{p) and Eq(jp) the dispersion relations of the fermions and gluons on the lattice 



We note that, in physical units, the potential V(Jl) (or its Fourier components V r ) are multi- 
plied by l/(a t al) = l/(^ s af). To avoid confusion we shall denote the (dimensionful) physical 
potential and its Fourier components by a superscript 'phys'. 

B. The meaning of r, the embedding of spacetime into color space, and why should 
there be a cascade of transitions with increasing fermion mass. 

Eguchi-Kawai spacetime reduction can be justified in perturbation theory. To see this 
one performs weak coupling calculations in the field theory and in the reduced model, and 
shows that the planar diagrams of the two theories are the same Since the fields in the 
volume-reduced theory have no spacetime coordinates, and consequently no momentum, one 
may wonder how can such an equivalence be valid. Put differently, what is it that flows along 
the lines in the Feynman diagrams of the reduced theory? 

The answer is tied with the fact that, in the reduced model, perturbation theory is per- 
formed around a particular ground state which is invariant. This means that the zero 
modes — the classical values of the gauge fields (or the logarithm of the Polyakov loops) 
- have eigenvalues that are uniformly distributed along the unit circle. If we denote the 
eigenvalues of these matrices by 9^ (here a is a color index and \x the index of the euclidean 
direction - we consider the generic case where more than one direction is compactified) , then 
the propagator of the fluctuating adjoint fields in the reduced model (the gluons and the 
adjoint fermions) would depend on 9 a . This dependence can be shown to be a very special 
one. The combination {9 a — 9S) behaves like the fi th momentum component of the (ab) th 
color component of the adjoint field. This is anticipated since in a gauge theory a momentum 
variable is replaced by a covariant momentum, and, in the presence of a classical background 
of the form A a ^ oc 6^9^, the covariant momentum of an adjoint particle is indeed (0® — 9^). 




\ 



J 



(2.7) 



(2.8) 
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If the ground state is invariant, then, in the large- N limit, the latter difference of angles 
would be uniformly distributed between —it and +7r, exactly what we expect from a lattice 
momentum variable. 

The fact that (0" — 9 ) behaves like a momentum variable can also be seen from Eq. f 12 . 5 j) . 
There, this difference is replaced by the variable u> G [— ir, 7r), and we see that the inverse 
gluonic and fermionic propagators (the arguments of the logarithms) are precisely the propa- 
gators of fields defined on R 4 and carrying a euclidean four momentum equal to (p, u). Since 
the winding number r is the conjugate to u (again see Eq. (12.51) ). then r behaves like a eu- 
clidean distance along an uncompactified dimension which replaces the single site that was the 
S 1 . 

The fact the winding of the Polyakov loops takes the role of a euclidean distance is most 
clearly seen from Eq. (12. 6j) . There r multiplies the energies of the gluons and fermions. In 
fact, in that equation, the different contributions to V r look like correlation functions of gluons 
and fermions that are separated distance r in euclidean space. 

Once we understand that the coefficients V r are determined by the difference between 
fermion and gluon correlators along a euclidean distance r, we can anticipate the cascade of 
transitions mentioned above. To see this note that when the fermions have a nonzero mass 
and are heavier than the gluons, they propagate a shorter distance - their contribution at 
large enough r is negligible compared to that of the gluons and only the latter will determine 
the behavior of highly wound loops. Thus, as we increase the fermion mass, lower and lower 
winding of loops become determined solely by the gluons. Since we know gluons generically 
cause the winding Polyakov loops to condense (indeed observe the minus sign of the gluonic 
contribution in Eq. (12.6jl ). we anticipate a cascade of transitions, where lower and lower 
winding of loops condense with increasing fermion mass. 

III. THE MINIMUM OF THE ONE-LOOP POTENTIAL 

In this section we minimize the one-loop potential. In Section IIII Al and Section IIIIBI we 
prove analytically that at m = and for < £ < 2 the ground state is invariant. Then, 
in Section fill CI we numerically calculate the values of V r for nonzero mass, for Nf = 1/2, 1, 2, 
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and for values of £ in the range [0, 4]. 



A. The massless case on the lattice 

In this section we ask what is the sign of V r when m = 0. We do so by comparing Ep 
and Eq- From Eq. (12. 6 p we see that V r will be positive for all r > 1 if Nf > 1/2 and 
Eg{p) > Ep(p). From Eqs. (I2.7l) - (l2.8p we see that the latter inequality will be fulfilled when 

P > 7— % 2 • 3.1 

1+ 2|P 2 



- - P • 7/' 

l2 



Using the identity p = p 2 — 4 X^ =1 sm4 (Pi/2), this inequality becomes 

^ 4 (|-l) >-4Esin 4 fe/2), (3.2) 

which holds if < £ < 2. Therefore, in this regime one finds that for Nf > | 

K > 0. (3.3) 

In particular, for one or more flavors we get V r > which means that V(f2) is minimized when 
tr [Q r ] = for |r| > 1, and the ground state is Zn symmetric. 

B. The massless case in the spatial continuum limit of a s — > and at fixed. 

Let us check what happens in the spatial continuum limit, where we keep a t fixed, but send 
a s to zero. In Appendix [B] we show how to take the £ — > limit of Eq. (12. 6ft . We find that the 
values of V r for r = 1, 2, 3 have quadratic, linear, and logarithmic UV sensitivities, respectively. 
The origin of such divergences was explained in Ref. 10] and calls for the addition of counter 



terms to the one-loop potential that are of the form |tr Q r \ 2 ; r = 1, 2, 3. These counter terms 
can always be chosen such that the renormalized V r= i^,3 obey 

Vi, 2 ,3 > 0. (3.4) 

Indeed, from the general form of Eq. (12. 6 j) with Ep(p) = Ec{p) = E(p), we see that for Nf > 1 
the choice in Eq. (13 .4p is a most natural renormalization condition. The rest of the V^s are 
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finite and in the massless case they have the following form. 

K P >f ■ = ^± f OrY z~ rm ; ^(P) = 2sinh- 1 b1. (3.5) 

(Note that here the integration variable p is dimensionless) . In Appendix [B] we show that 
Eq. (13. 5p is what one obtains if one chooses to use a continuum regulator for the spatial UV 
(for example, a hard cutoff or dimensional regularization). Interestingly, while the integral in 
Eq. (13. 5p is the correct £ — )• limit of Eq. (I2.6P only for r > 4 (otherwise it is UV divergent), 
then we note that inspecting the UV behavior of the right hand side in Eq. (13.51) shows that it 
is finite for all r > 2, and UV linear divergent only for r — 1. This linear divergence is one of 
the UV sensitivities of the one-loop potential that Ref. exposed by analyzing a continuum 
regulator. 4 This means that if we regulate the spatial UV of the single site theory on a lattice, 
and then take the £ — > limit, then the UV sensitivity of V phys '(Q) is more severe than what 
we find in a continuum regularization of the spatial UV. For further details on how to obtain 
Eq. (13.51) from continuum regulators see Appendix [Ql 

The simple form of the dispersion relation in Eq. (I3.5P allows us to analytically calculate 
the coefficients V phys \ We do so by inverting the dispersions, \p\ = 2sinh(.E/2), and writing 
Eq. (E31) as 

K p hy s. = ( 2N f - 1) r dE sinh (£/ 2 ) S i n h(£) e~ rE . (3.6) 
n 2 raf Jo 

For r > 2, Eq. (EID leads to 

T/phys. _ 2N f ~ 1 1 (a 7) 

(r 2 -i)(r 2 -!)' { } 

which is positive for Nf > 1. Combining Eq. (13.41) and Eq. (13. 7p we see that the one-loop 
potential for Nf > 1 and m = is minimal when the ground state respects the symmetry. 

As a check of our calculation we note in passing that for large r the right hand side of 
Eq. (13. 7p is precisely the one that appears in the calculation of Kovtun-Unsal-Yaffe in Ref. {(| 



4 Another divergence was argued to be present in the r = 2 term, but here we see that when we regulate the 
compact dimension with Wilson fermions that have a 'Wilson parameter' equal to unity, the coefficient of 
the linear divergence is exactly zero. 
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(if we identify the distance a t with the extent of the compact direction in that calculation). 
This is expected from Eq. ( 12. 6p . There, we see that when r is large, the spatial momenta 
that would contribute to the integral have a corresponding energy that obeys E < 1/r -C 1. 
For such momenta we can replace E by 2sinh(i?/2) = \p\, and we arrive at the continuum 
equations that appear in Hef/a. T^ looW at la Ig e^n gr >> 1 leads — ytoth e 
continuum limit in the compact direction. We extend this observation to the massive case in 
Appendix [Cj where we show that in the large-r limit V^, phys ' agrees with the result of Ref. 8j. 

In Appendix [D] we perform an additional check of our calculation and verify that Eq. (13. 7p 
is obtained with a different mathematical procedure (i.e. first integrating over p and then over 
uj) and within a different regulator for the spatial R 3 directions. Specifically, we compare with 
a hard cutoff and with dimensional regularization. 

Finally, we note that the continuum limit of our results for massless quarks is also consistent 
with the single-site case of Ref. 17J (There the non-compact directions were regularized in 
the continuum, and the compact direction was allowed to consist of T > 1 sites). 



C. Minimizing the lattice one-loop potential for general values of m, £, and Nf. 

We calculated V r for general values of a t m, £, and Nf by numerically performing the spatial 
integrals in Eq. ( I2.6p . 5 For each choice of parameters we denote by K the smallest value of r 
for which V r is negative. For this parameter choice breaks to Z%. Since we numerically 
calculated V r for r = 1, 2, 3, . . . , 20, then, from these results, we can only conclude about 
symmetry breaking to subgroups that are at most Z20. Nevertheless, in Appendix O wc 
investigate the large-r limit and show that in that limit the symmetry realization is identical 



to the one of the continuum calculation done in Ref. 



In particular, we see that the when 



a t m <C 1, there is an instability for Z N to break to Z K with K 1 if 

a t m>c Nf /K ; c lj2 ^ 2.03, 3.155. (3.8) 



5 The numerical integrals were performed with a trapezoid method and a resolution of the interval [— 7r, 7r) in 
each direction equal to 2ir/L with L = 30, 50, 90, 130, and 150. The results we present were stable under a 
change in L. 
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For generic values of a t m, where the breaking of is to with K ~ 0(1), the values of 
the lattice V r 's (and the resulting phase diagram) can be obtained by numerical integrations 
of Eq. (12. 6p . and we present our results in Figs. ([THS]). 




FIG. 1: Large- N phase diagram for theories with a single Majorana Dirac fermion. 

As the figures show, for all values of £ we find that the ground state at m = has the 
highest possible symmetry. Since we restricted our numerical study to instabilities of Zn 
breaking to Z r with r < 20, we cannot conclude whether Zn breaks to Z r with r > 20, but 
from Eq. (13.81) we see that at large enough values of r and at a t m = 0, then the Z^ is intact 
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FIG. 2: Large- N phase diagram for theories with two Major ana fermions. 

for all £. 

When a{m increases from zero, the Zn symmetry breaks to a subgroup Z r with r a de- 
creasing function of a t m. We emphasize that while the figures seem to indicate that there 
is a finite range of a t m for which the in intact (see white patch just above the x-axis in 
Figs. fl2H3])), then this is most likely an artifact of our restriction to investigate numerically 
symmetry breakdowns of the form Zn — > Z r with r < 20. In particular, as we mention above 
(see Eq. (13. 8ft ). in Appendix O we see that even for infinitesimal values of a t m, there is an 
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FIG. 3: Large- iV phase diagram for theories with four Majorana fermions. 



instability for the symmetry to break to a subgroup Z r with r ~ (a t m) 1 . 



IV. SUMMARY 



The calculation in this paper is concerned with large- N lattice QCD with adjoint fermions 
on R 3 x S 1 . We calculate the one-loop potential V of this theory as a function of the Polyakov 
loop Q that wraps the compact direction. We use Wilson fermions and the compact direction 
consists of a single site. Our calculation extends Ref. [lOj by that it directly determines, for 
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different values of the quark mass, the lattice anisotropy parameter £, and the number of 
flavors, Nf, the instabilities in V(Q) toward the breakdown of the Zn center symmetry to a 
Zk subgroup with K < N. 

Our results are summarized below. 

• For m = we prove analytically that the minimum of the one-loop potential respects 
the Zn symmetry. This is true for Nf > 1 and in the regime < £ < 2, which includes 
the case where one takes the continuum limit in the uncompactified R 3 directions. 

• For m = and Nf = 1/2 (a single Majorana fermion) the one loop potential is zero 
in the continuum limit of the uncompactified directions and so one cannot use it to 
determine the ground state. For any finite a s the one loop potential is nonzero and 
prefers a Zn symmetric ground state. 6 

• For m > we see a cascade of transitions where the Z^ breaks to a smaller Zk subgroup 
with K < N, and K being a decreasing function of a t m. At a t m = oo the center Z N 
completely breaks down. 

• When we investigate the possibility to break Z^ to Zk with K ^> 1 with the lattice 
one-loop potential, we find the continuum results of Ref. jg|: when m is increased to 
values of 0(1/K), then Zn breaks to Zk- Thus the Polyakov loops with large winding 
numbers behave like they do in the continuum. 

The fact that we see symmetry breakdown of Zn for a wide range of masses means that one 



of the conclusions of Ref. ICj does not hold. There, we concluded that the Zn symmetry is 
intact even for moderately heavy masses, but relied on comparing the one-loop potential only 
for ground states that either leave Zj^ intact or break it down to or Z 2 . Here we see that 
for intermediate masses the symmetry breaks to Z r subgroups. 

We wish to emphasize that, while our results are very similar to those obtained in the 



continuum works Refs. [6|-|8j, then the UV sensitivities of the one- loop potential discussed in 



For a discussion on the nonperturbative effects of instantons see Ref. 
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Ref. 



lOj make it impossible to determine a priory whether this similarity will take place or 



not (this is because our lattice theory has only one site in the compactified direction). 

The methods we use make clear that the winding number of the Polyakov loop around the 
compact direction is identified with a euclidean distance along the direction one reduces (at 



least in the single-site case that we studied in this paper). As we explain in Section III Bi this 
leads one to anticipate the cascade of transition seen in our calculation and in Refs. (7J, |8j: 
when m > the fermions are heavier than the gluons and the winding- number /distance 
identification tells us that their effect will propagate to smaller number of windings than that 
of the gluons. Therefore, at nonzero m, the behavior of large enough winding numbers will be 
determined by gluons, and these will cause a condensation of the multiply wound Polyakov 
loops, leading to a partial breakdown of the center. 

Finally, we wish to remark on the complications of inferring from the results obtained in this 
paper on nonperturbative lattice simulations. First, we work at weak coupling, and it is not 
clear how the stronger fluctuations induced by stronger lattice couplings (where one would 
actually perform lattice simulations) change the one-loop picture. Second, our spacetime 
geometry is that of R 3 x S 1 , while standard lattice simulations are often done on a four-torus. 
Indeed, the one-loop potential for a system that has more than one reduced direction depends 
on more than one Polyakov loop. Furthermore, if one defines the theory on a lattice with 
a single site in all directions, then the one-loop calculation will not be valid; there are zero 
modes around specific vacua ('singular tolerons' - see Ref. [uj]) that cause IR divergences and 
require a nonperturbative treatment. 

The last difficulty that arises when one tries to infer from our results on what would happen 
in simulations is related to 1/N corrections. Here we note that when we minimized the 
one- loop potential we assumed that all the windings of the Polyakov loops are independent; 
an assumption that holds only when N = oo. At finite values of N the minimization is 
more complicated and, for the continuum one-loop potential, was performed numerically in 
Refs. (3, S|. The result of these references is interesting: if we denote the size of the S 1 by 
L, then one finds that there is a critical value of (m x L) = fi c , below which the cascade 
of transitions ends and the Zn symmetry is completely intact. In particular, just above /i c 
the Zn breaks to its largest subgroup, Z N / 2 . Thinking about this breakdown as Zn — > Zk 
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with K ^> 1, we can anticipate that fi c {N) ~ 1/N (see last bullet above). Thus, while for 
given values of N, L, and m that obey \l < fi c , the might be intact, then, according to the 
R 3 x S 1 one-loop analysis, an increase in N (done at fixed m and L) will result in a breakdown 
of Zpj. Since lattice simulations are always done at finite value of N, then the breakdown we 
discuss above and the cascade of transitions we discuss in this paper, might occur only if one 
systematically checks how the simulation results change when N is increased. Such systematic 
checks were done in Ref. [l^] to the extent possible with the metropolis simulation algorithm 
used in that paper (in practice this algorithm allowed the study of 8 < iV < 15). It will be 
very useful to develop faster algorithms and simulate larger value of N, so to extend these 
systematic checks further. 

The caveats mentioned above might explain why the nonperturbative lattice studies pre- 

nnn 

sented in Ref. [13|, H5J, H6( observed a Zn symmetric ground state even at nonzero, but small, 
mass. 7 8 
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Note added 



For m — 0, the recent preprint 17| by Poppitz and Unsal 
to the one we study in this paper. Specifically, Ref. [12] 



is studying a similar construction 
defined the theory to lie in R D x S 1 , 



15] finds an elaborate phase structure for SU(3), it reports a center invariant phase at sufficiently 



16j used staggered and overlap fermions, respectively, and so the results we obtained in 



7 While Ref. 
weak couplings . 

8 Ref. and Ref. 

this paper do not directly related to that work (the one-loop potentials on lattices with a small number 
of points sometimes have UV sensitivities that requires the addition of certain relevant operators when 
comparing results obtained with different lattice actions.) 
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and regularized the UV of the R D in the continuum and the UV of the S l on a lattice of V sites. 
For the case D = 3 and T = 1, this construction becomes the one we study in Section fill Bl 
- i.e. the limit in which we send the lattice spacing in the non-compact directions to zero. 
Indeed, in that limit, our formulas for V r (see e.g. Eq. (13.71) and those in appendix [D]) become 
consistent with those of Ref. 1171 1. 



Appendix A: The calculation of V r . 



In this appendix we show how to arrive at Eq. (12. 61) . We begin by simplifying the argument 
of the fermionic logarithm of Eq. (12.51) and write 

. 2 



■^p + sin 2 u + m 2 v (p,u) = 4sin 2 u;/2 (1 + a t m + T^Pj + + ( a t m + w?P 2 ) > ( A1 ) 

1 „o\ . . o JzP 



2f 



1 " 2 / 1 -9> 2 

p p + [a t m + Yf p A 



1 + a t m + ^-f I 4 sin 2 u/2 + 

2 i / ! (l + a t m+±p 2 ) 



(A2) 



We factored out ^1 + a t m + \p 2 ^j because it is independent of u and does not contribute to 
V r when r > 1. We can now write 

4 sin 2 (u/2) + 4 sinh 2 (E G (p) /2) 



-2N f log [4 sin 2 (u/2) + 4 sinh 2 (E F (p) /2)] } , for r > 1. 



(A3) 



with Ef : g(p) the lattice dispersion relations of the fermions and gluons that we give in Eq. (12. 7p 
and Eq. (|Z%jl . Defining I r [E] as 



we have 



/,.[£'] = Re / ^ e irw log [4 sin 2 u/2 + 4 sinh 2 E/2 



V r = J {^j {I r [E G (p)}-2N f I r [E F (p)]}. 



(A4) 



(A5) 



To proceed, we need to calculate I r [-E]. We do so by first integrating by parts, then performing 
the change of variables z = e tul , and finally using the Cauchy integral theorem: 



I r (E) = Re 



d (e iruJ \ 
2-ftir 



log 4 sinh 2 E/2 + 4 sin 2 u/2 
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Re log 4 sinh 2 E/2 + 4 sin 2 u/2 - / 



2vr /- 27r gL> , 2sina; 



2mr 4 sinh 2 E/2 + 4 sin 2 w/2 



/ 1 / dz z r+1 — z r ~ l \ e~ rE 
= Re {-r?2^ ( Z -e-z)(z-eZ) ) = - — - (A6) 

Combining Eq. ( 1A6j) and Eq. (1A5[) gives Eq. (12. 6p . 



Appendix B: The £ limit of Eq. 

In this technical appendix we show how to take the £ — > limit of Eq. (12. 6ft . To do so we 
divide the integration over the Brillouin Zone into two regimes. 

Regime I. The neighborhood of the Brillouin Zone origin. Here p ~ p ~ |p| = £ g with 
g ~ 0(1). It is straight-forward to show that the contribution of regime (I) to 
the one-loop potential is given by Eq. (I3.5P : 

VST = f (v) 3 e ~ rm ! E(p) = 2mhbT 1 \Pi. (Bl) 

H a t J-oo \27T J 

Regime II. The rest of the Brillouin Zone, where p ~ 0(1). In that regime we have 

E G m 21og(p/0, (B2) 



E F ^\o g [[ 2 f + \p*\/i\. (B3) 



Substituting Eqs. (jC2tiC3]l into Eq. (J2l| (and dividing V r by £ 3 c4 to get 1/ phys ) 
we find 

F p hys . regime (II) l_ ^ Q ^ + ^ ^ 2r _ 3 ^ ^ 

Where here the first term is the contribution of the fermions, and the second of 
the gluons. 

Combining Eqs. (IB II) and (IB4I) we find that 

• The gluons render the V r =\ linearly divergent. 

• The fermions render K=i,2,3 quadratically, linearly, and logarithmically, divergent, re- 
spectively. 

• Only V r >4 are finite in the £ — > and are given by Eq. (13.51) . 



Appendix C: Comparison of l/ phys in the r — > oo limit with the result of Ref. 

In this appendix we show how, for large values of r, the coefficients V r in the lattice one-loop 
potential, become those obtained previously in the continuum calculation of Ref. [8j. This 
serves as a check of our results and also allows us to conclude that for small, but finite, values 
of m, there is an instability for the to break to Z% with K ~ 1/m ^> 1. 

Our starting point is to take the large-r limit of Eq. ( 12. 6ft . For r ^> 1, the energies Ef,g{p) 
that contribute to V r are those that obey 

E f ,g(p) % Vr < !■ (CI) 

For such small energies we can write 

Eo(p) = \P\ (C2) 
MP) = {\p\ 2 + (Ca t m) 2 ) 1/2 /C. (C3) 

Substituting Eqs. (1C2HC3I) into Eq. (12.61) and using straight-forward manipulations gives 

V^- = - i } T1 (N f tf.KM-i) ; fi r = ra t m. (C4) 



Identifying a t with the length L of the compactified direction in Ref. [6j and Ref. js], we see 
that the large-r result in Eq. (jC4p becomes identical with the results of these works. 

The meaning of Eq. (12.61) is that when the expression in the brackets becomes negative, 
there is an instability of the symmetry to break to Z r . This will happen when jj, r obeys 

fi r > 2.03 N f = 1 (C5) 
fir Z 3.155 N f = 2. (C6) 

This gives Eq. (13. 8p . 



Appendix D: Comparison of the result obtained in the continuum limit of the un- 
compactified directions (Eq. (13. 7|) ) with other regulators of the spatial UV. 

As a check we compare the result obtained in Eq. (13.71) with the one obtained employing a 
continuum regulator for the spatial / d 3 p integrals. The regulator we choose is a hard cutoff 
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in momentum space, but we also remark below on what one would obtain with minimally 
subtracted dimensional regularization. 

When the spatial space is a continuum, the one-loop potential of the massless theory is 
given by (here the factor of l/a t ) is required to make the potential dimension four) 



V(Q) 



a t 



4 2 (d a -6 b ^ 
-tt sin 



a 



where I [M 2 ] is 



'dk\\ 



k 2 + M 2 



(Dl) 



(D2) 



This integral has both cubic and linear divergences. The cubic ones can be removed by 
subtracting 1(0), and with a hard cutoff we obtain 



I(M 2 ) - 1(0) 



6tt 2 



'dk\\ 



d(k 3 ) log 



M 2 
k 2 



1 



2tt 2 



k 2 dk log 



M 2 
1 + ^ 



1 + 



M 2 



6tt 2 



A; 3 log 



A; 3 M 2 (-2) AT 
1 + ^ 



1 + 



M 2 



AM 2 M 2 



AM 2 |M| 3 r A /\ M \ dx AM 2 \M\ 3 f A/m dx 
6ir 2 3% 2 Jo 1 + 1/x 2 6n 2 3n 2 Jo 
\M\ 3 rNW\ dx AM 2 AM 2 \M\ 3 _ x 



3tt 2 Jo 1 + x 2 6ir 2 



+ 



3vr 2 



3tt 2 



tan _1 (A/M) 



AM 2 \M\ C 



(D3) 



2vr 2 6tt 

In the case of dimensional regularization (which implicitly sets the linear divergence to zero) 
we get 



I(M 



2\ dim— reg. 1-^1 



67T 



(D4) 



Let us now calculate the moments V r of V(Q). To do so we start from Eq. (ID3|) which we 
write as (here we drop constants that do not depend on 9) 



V(Q) 



'1 - 2N 



fJ 



E 



;i - gjy>) 

2?ra t 4 



2a t A 

7T 

2a 4 A 

7T 



\ 8 


(e a -e b \ 


3" 


) + 3 


Sm [ 2 J 











trfi| 2 + ^ 



a^b 



'e a - e b 



sin 



(D5) 
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If we write 



then we have 



|sin (w/2)| 3 = 2 W r cos (ur) + W , 



V r 



2N 



f 



2naf 



r=l 



2*A + 

7T ^3 



for r = 1. 
for r > 2. 



(D6) 



(D7) 



The calculation of W r involves a standard Fourier transform: 

™ sin 3 W /2 = -— Re / e l ™ * 

7 8tt y 

1 .. -i(-l) r -l i(-l) r -l 



Re 



Re 



87T 
1 

'8tt 



-i(-iy - 1 r - i 



i 2 (r + 3/2) 
-3 3 



which gives 



9/4 



+ 



1/4 



i 2 (r + 1/2) 
3 

+— 



z 2 (r-l/2) « 2 (r-3/2) 
1 



4tt (r 2 - f)(r 2 



(2iV/ - 1) 
7T 2 af 



at A 
2 



16 
15 ' 



for 



1. 



for r > 2. 



(D8) 



(D9) 



(r 2_| ){r 2_l) 

Eq. flD9|) is precisely what we obtained in Eq. (ID 1 p for r > 2. 

The fact that V\ is UV sensitive is explained in Ref. [10] and means that the naive continuum 
limit one takes to arrive at Eq. (1D1I) requires a counter term. This counter term can be chosen 
such that the renormalized V\ has arbitrary sign, and so such choice can affect the realization 
of the Zjy symmetry. It also indicates that the realization of the lattice theory can depend 
on the particular action one chooses to work with. For further discussion on this point see 
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